Ci 14 
(310/1961) UDC 539.3:624.074. 4:69.024.4 


ACTA 
POLYTECHNICA 


CIVIL ENGINEERING AND BUILDING CONSTRUCTION SERIES No. 14 


OSMUND JENSSEN 


Shallow Hyperbolic Paraboloidal Shells 


Norwegian Contribution No. 14 
Trondheim 1961 


> 


ACTA POLYTECHNICA SCANDINAVICA 


. @ Scandinavian contribution to international engineering sciences 


Published under the auspices of the Scandinavian Council for Applied 
Research 


in Denmark by the Danish Academy of Technical Sciences 


in Finland by the Finnish Academy of Technical Sciences, the Swe- 
dish Academy of Engineering Sciences in Finland, and the State 
Institute for Technical Research 


Norway by the Norwegian Academy of Technical Science and the 
Royal Norwegian Council for Scientific and Industrial Research 


Sweden by the Royal Swedish Academy of Engineering Sciences, 
the Swedish Natural Science Research Council, and the Swedish Tech- 
nical Research Council 


Acta Polytechnica Scandinavica is divided into the following sub-series: 


Chemistry including Metallurgy Series, Ch 

Civil Engineering and Building Construction Series, Ci 
Electrical Engineering Series, El 

Mathematics and Computing Machinery Series, Ma 
Mechanical Engineering Series, Me 

Physics including Nucleonics Series, Ph 


For subcription to the complete series or to one or more of the sub-series and 
for purchase of single copies, please write to 


ACTA POLYTECHNICA SCANDINAVICA PUBLISHING OFFICE 
Box 5073 

Stockholm 5 

Sweden Phone 67 09 10 


This issue is published by 
NORGES TEKNISKE VITENSKAPSAKADEMI 
Trondheim, Norway 


and is designated as NTVA Publication series 1, No. 12, 1961. 
Editor: Knut Alming 


7 
| 
H 
3 


Ci 14 
(310/1961) UDC 539. 3:624.074.4:69.024.4 


ACTA 
POLYTECHNICA 


CIVIL ENGINEERING AND BUILDING CONSTRUCTION SERIES No. 14 


OSMUND JENSSEN 


Shallow Hyperbolic Paraboloidal Shells 


Norwegian Contribution No. 14 
Trondheim 1961 


| 
4. 
r 


Aktietrykkeriet i Trondhjem 


i 


SHALLOW HYPERBOLIC PARABOLOIDAL SHELLS 


BY 
OSMUND JENSSEN *) 


SUMMARY 


This paper deals with the bending theory of shallow hyperbolic paraboloid shell 
roofs. The solution is obtained in terms of power series, which converges for extremely 
shallow shells. The method to obtain the solution in the divergent domain is indicated 
as well, 


Fig. 7. 


INTRODUCTION 

In this paper the bending theory of Marguerre [1] for shallow shells is applied to 
the shell roof shown in fig. 1. The two differential equations of this theory are reduced 
to one single complex equation. This equation is separable but the Levy-type solution 
is not too well suited for the involved boundary conditions of our problem. In this 
work the differential equation is transformed into an integral equation in the way 
shown by Vekua [2]. The solution is obtained in terms of power series. Further the 
boundary conditions are established and an iteration procedure to satisfy these condi- 
tions is discussed. 


BASIC EQUATIONS 
The unstrained middle surface of the shell in the first octant fig. 1 obeys the equation 
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For this surface Marguerres [1] equations for the displacement w in z direction and the 
stress function @ are given by 


D Véw + 2k Indy = ds (3) 
1 Yew 
— 2k —— =0 (4) 
where 
Eh3 
C = Eh, 


and q, is the load intensity in z-direction. C is the stiffness in tension multiplied by 
(1 — v®) and D the stiffness in bending. E is Young’s modulus, v is Poisson’s ratio 
and h is the thickness of the shell. 


Multiplying eq (4) by an arbitrary constant ko and nay it to eq (3) give 


+ ko] — 2k koC [w— che — (6) 
The expressions in the brackets will be equal if 
1 
giving 
— 
Eh? 
(DC) 
With this value of ko equation (6) is transformed into 
YU 
ay — on 
=D (8) 
where 
U=w+— (9) 
(DC) 
and 
e = ik (10) 


Forces and moments are expressed by w , ® and the load intensities qx , qy in x and y 


directions as follows 
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These forces and moments shown in fig. 
(2a-c) are applied to an element cut out from 
the shell by the planes x = const. x + dx = 
const., y = const. and y + dy = constant. 
Nx, Nxy and Ny are tangential to the midd- 
le surface, hence their vertical components 
are included in the effective vertical edge for- 
ces Rx and Ry. 

The strains of the deformed middle surface 
are 


1 
= rol [Nx — 


ey = [Ny— NG] (12) 


As shown by W. Fliigge and F. T. Geyling [3] these strains may also be expressed in 
terms of the first derivatives of the displacement components u, v and w in x, y, z 
directions respectively as follows 


from which we obtain 


& = + 
wv 
dy 
du 
Y= dy 
by 
dy? by 
_ 
Ox 


(13b) 


These equations are needed to establish our boundary conditions. 
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INTEGRATION OF THE HOMOGENEOUS DIFFERENTIAL EQUATION 


In the equation to be solved, namely 


YUn 
4 
V Un (14) 


where index h indicates that we deal with the homogeneous equation, we introduce 
the new complex and dimensionless variables 


in (15) 
C= 
where 
fa (16) 
and obtain the equation 


By integration this equation takes the form 


Un=—* Uf — SS (18) 


as easily obtained from eq (17) by use of eq (2) (10) and (5). For Un in the right side 
of the eq (18) we now substitute the biharmonic polynoms 


with 


Unn = + + + (20) 


where An, Bn, Cn and Dy are arbitrary constants. 
Let U! be the function obtained in this way on the left side of eq. (18). Next we 


substitute U! | on the right side of eq (18) and obtain on the left side the new function 
Uj, : Going on in this way we find 


q=0 


—1)! 
2q)! Bop 1) Haas} (22) 


‘ 


where 


Page + (— (23) 
and 


a= + (— 41 


The polynoms U¢) in eq (22) should satisfy the equation 


ane +* er] = ° on 
which can be verified by direct substitution. 

It is often convenient to work with solution wich are symmetric or antisymmetric 


with respect to € and 7. These are obtained by proper choice of the constants of 
integration as given in the following table. 


n Symmetric solution Antisymmetric solution 
n = 4j A B : Cc D i! c A B ‘2 b Cc D ; d 
= _ = = i- = i- 
4] a a 2 a n 2 a 2 2 
45 +2 A B Cc D A B Ca=D i! d 
n= 1 = = i- a= = — = a 
A] a n > n 
1 +i 1+i i— 
Aa= 2 an = Ch As = ba Ca = = 
n=4+1 
1 1—i 1— 1 
1+i 1—i 1— 1 
n= 4) +3 
1+i 1 


where n, and j are integers and an — dn are new complex constants of integration. 
The solution given by eq (22) and (23) is independent of the choice of origin, hence 
we may add the solution obtained by putting the origin at the point & = 1 and = 1, 
which corresponds to the substitution of § — 1 and 7 —1 for & and » in eq (22) and 


(23). Let the constants of integration connected with this solution be denoted by 
a—d.. 
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INTEGRATION OF THE NON-HOMOGENEOUS EQUATION 
The considered equation namely eq (8) may by use of eq (16), (10) and (19) be written 


ata, 


When U is a smooth function and x a large parameter, we can neglect the fourth deri- 
vatives in this equation and obtain 


YU, a‘ 


where index p indicates that we deal with the particular solution. For q,; = qo = const. 
we obtain from eq (26), (9), (11) and (19) the following well known membrane solution 


YO, 1 ab 


= = — at (27) 


In this case and when q, is a linear function of & and 7 eq (25) is identically satisfied. 
For other loading functions and smaller values of x the solution can be found by the 
method of successive approximations. 


BOUNDARY CONDITIONS 


We consider the shell shown in fig. 1 with qx = qy = 0 and qz = qo = const. 
applied to the entire roof. The edge stiffening members for x = a and y = b are par 
example assumed to be rigid in the direction of their axes and have negligible bending 
resistance in planes tangent to the shell. Hence w , the edge moment, Nx and Ny should 
vanish along the outside boundary giving 


which can be written as 


U= VU =0 (28) 


For y = 0 we have from symmetry the conditions 


(29) 
and y 
y= 0 together with yx_y_o = 0 


which according to eq (13), (12), (1) and (2) can be written 


| 
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beat 
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dev 1 (2(1 + ») INxy oN, k 


and 


(Nxy)x-y.0 = 0 (30b) 


This edge has to be stiffened by an edge beam with area and moment of inertia 
denoted F, and I,. Hence we have to add the two conditions 


(€x)y.0 = &0 (31) 
and 
32) 


where the terms with index zero are the strain and curvature of the edge beam, We 
now assume that the neutral planes of the edge beam and the shell coinside at the 
line y = » = 0. Hence the Nxy forces on both sides of the edge beam give 


and equation (31) takes by use of eq. (12) the form 


a 
(ex)ya0 = [Nx— =, (33) 


The vertical load pr. unit length of the edge beam consists of its own weight and 
external load denoted by g, and the effective edge force Ry on both sides of the edge 
beam as shown on fig. 3. This gives the shearing force 


| ~ Q=— J (2Ry + g,)dr (34) 
0 


N, 
and the bending moment 


om 


fig. 3. ade =f de J + (35) 


which determines the curvature of the beam given by the equation 


d2wo M 
(36) 
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Eq (32) may then by use of eq (34-36) be written 


a 
EL de J (2Ry + g,)dx = 0 (37) 


Our boundary conditions at the edge y = 7 = 0 given by eq (29), (30a), (33), and 
(37) can now be given in terms of U by use of the equations (11), (9), (1), (16) and (19) 
as follows 


R (5) 0 
min do fe + - 0 


where Re [ ] and Im [ ] denotes the real and imaginary part of the term in the brackets 
respectively and 


EI 2Ca 


Da’ % EF; 


a; = 


The similar expressions for the edge & = 0 are obtained by substitution of 7 for & 
and & for 7 in eq (38) and (39). 


These equations have to be supplemented with eq (30b) which can be written 


Im = 0 (40) 


The boundary conditions are expressed by the aerivatives of U with respect to & and 7 
while the homogeneous solution is expressed in terms of ¢ and £, hence we list the 
following transforms where m = n + 2q — 4a 
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(n 4+ 2q — 2a) 4- (— 1)9ByEm="] 
(2a) [Ang™ 4 (— 
=i 
i {(n 4+ 2q — 2a) 4 (— 


(n 4 2q — 2a)(n 4 2q — 2a — 1) +4 (— 
4 2(2a)(n 4- 2q — 2a) [Ang™ + (— 1)9BnE™] 
+ (2a) (2a — 1) [Ang™*? (— 1)9ByEm +2] 


= [aga + 2 — Powe 


—(n + 2q — 2a) (n + 2q — 2a — 1) + (— 
2(2a)(n 2q — 2a) [Ant™ + (— 1)9Bn€™] 
— (2a)(2a — 1) [Ant™*? (— 
= i {(n + 2q — 2a)(n + 2q — 2a — 1) +-( — Bal m — 2 
— (2a)(2a + (— 1)9+1B,Em+2] 
((2 + ») =i [a + ») + (5 + ») —(5+ 
—(1 + ») Pose 


which in turn is eqnal to the expression on the following page 
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i (1 4+ »\(n + 2q — + 2q — 2a — 1)(n 4+ 2q — 2a — 2) 


(5 + 2q — 2a)(n 2q — 2a — 1) 4 (— 
(5 v\(2a)(2a — 1)(n 4+ 2q — 2a) (— 
~(1 + »)(2a)(2a — 12a — 2) +4 (— 1)4 


The similar expressions for the polynoms Hyqs are obtained by substitution of 
Hnap for Page ’ Cs and Dy, for An and Ba and 2p | 1 for 2a. 


APPLICATION OF THE SOLUTION 


We take as an example a shell witha = b, a, = a, 4, B, 
v = 0 and with uniformly distributed load q, — qo applied to the entire roof (fig. 1); 
hence we need only the symmetrical part of the solution eq. (22). 

We have to calculate two sets of constants namely ap, Cy and ap’ , Cy’ corresponding 
to the solutions U and U’ obtained when the origin of coordinates in eq (22) is placed 
at & = 7» = 0 and & = 7 = 1 respectively. These constants can be calculated from the 
equations obtained by equating to zero the coefficients of &" for each n at the line 7 = 0 
and of 7" where & = 1. This infinite set of equations will not be easy to solve. However, 
by neglecting the influence of the solution U’ at the boundary 7 = 0 we obtain a set 
of equations from which the constants ay and Cp are easily found. ‘This solution gives 
disturbances at the boundaries € = 1 and 9 = 1 which can be eliminated by addition 
of an U’ solution. This gives disturbances at the boundaries & = 0 and 7 = 0 which 
may be eliminated by a new U solution, which again has influence at the outer bound- 
aries, and so on. 


The convergence of this iteration procedure will depend on the Nx and Ny forces on 
the boundaries & = 1 and 7 = 1 obtained from the first U solution. The first U’ solution 
will then include equal and opposite Nx and Ny forces on these boundaries and when 
these forces have a resultant there will be a beam action between the parallel boundaries. 
This beam action is one of the main causes for the reflected disturbances on the inner 
boundaries. ‘The more evenly distributed these Nx and Ny forces is the more undamped 
will they go from one edge to the opposite with correspondingly bad convergence 
of this iteration procedure. Let us now consider the first step namely the calculation 
of the constants an , Cn from the conditions at the boundary 7 = 0 given by the equations 
(38). These conditions result in the following set of equations for the real and imaginary 
part of the coefficients denoted by subscripts r and i respectively 
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aor = 0 


3asr — Car = 0 


3 
5 1 
10asr — + 2 “agi = O (42) 
12a¢r +- = 0 


— + = 0 
36c4, + 28xcor = 0 
60a5; — + 48x(asr + Car) = 0 


120aq; + 48xagr — 


(43) 
+ y Im Gr), = 
— 3ag; + cai = 0 
— — + (3asi — Cai) = 0 (44) 


20as1 + + 3xasr — xCar + 7 (4ca + =0 


— 10ce: + — (Sasi — — + xca1) =0 


2acer — Re do f de + 
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agr + Car = O 


a(12a4r — — Baar + Scar) — = 0 


(49) 
a(20a5r + 20c5r — 3 nea) + (36C4r — 28x21) 


1 


a(30cer — ; wea) + — + — 0 


From these four sets of infinite equations all the constants with subscripts of higher 
number than two can be expressed by agi, Car and cz, and these constants can then 
be determined from the first of eqs (44) and (45) together with the boundary condition 
(40). By aid of (27) this can be written 


Im U a‘ 0 
dF dy l6xD 4° 
and by use of eq (22), (41) and (19) we obtain 
a‘ a‘ 


CONVERGENCE OF THE SERIES 


For large values of n and j the terms with the largest subscripts in the equations 
(42-45) will dominate and the constants may be calculated as follows. 

We want to calculate par example ayaj,2); . ‘This can be done by : me of the boundary 
conditions (38b) and (38d). By equating to zero the coefficient of &~' the first of the 
above mentioned conditions gives 


(4))® (1 +0 — +0 + (ar +0 
+ 8x(4j)* (1 +0 Gl — +0 Gl =0 


which has been developed by use of eq (22) and (41). When j is very large we may 


neglect the terms of order of magnitude a! and eq (47) results in 
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In eq (38d) the first term is proportional to (4))* while the next two terms are propor- 
tional to (4j) and 8 respectively, hence this boundary condition may for sufficiently 
large values of j be written as 


Re = 0 


(49) 


4j.2 


Equating to zero the coefficient of €~" in this condition we get by neglecting all 


terms of order of magnitude i the equation 
2 3 


x 


We now substitute this value of ac4jr in eq (48) and obtain 
1\* 1 | | 
In accordance with eq (48) we obtain 


2! 


which substituted in eq (51) gives 


1? 2 

In the same way eq (50) and (48) result in 

2 

or 4! 6! (54) 


For other coefficients we obtain similar expressions, where the main term on the right 
2 2 

side will be either (5) as in eq (53) and (54) or (5) . The series seem therefore to 

be convergent when x < 2 and divergent for x > 2. However, no singularity in the 

solution should be expected within or on the boundary of the shell, hence one of the 

methods of summing divergent series may be used, For example transformation into 

continued fractions after Ruthishausers «Q. D. algorithmus» [4]. 


To calculate cer and cg; we want the value of the function Im (5) where &= 1. 


The contribution to this function from the coefficients of &” +l can be written 
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j-0 p-0 

where the index 8 in ai4j,2-2p) 0 i8 equal to i, when p is an even number and equal 
to r when p is an odd number. ‘The series given by (55) will according to eq (53) and 
(54) for a given value of & result in a function with a singularity for a sufficiently large 
x-value. Hence, the contributions from the different types of coefficients can not be 
given separately and we must express all the different coefficients by external loads, 
that is by ag and f, before the series expansions of the different derivatives of the 
solution can be summed. This however is impossible without knowing the values 
of ca; and Cor. ‘These coefficients may for x > 2 be found in the following way. 


J 
In the first of the equations (44) we substitute for Im (\\) at the points = 1, 


» = 0 the contributions from the different types of coefficients. In these expressions 
which are similar to (55) we change the upper limit of summation from infinity to a 
finite number denoted by m. Further we treat the first of the equations (45) in the 
same way and denote the values of cer and cy; obtained from these equations (Cer)m 
and (ce:)m. ‘The values of cer and cz; may then be found by the equations: 


Cer ( 2 (C2r)m) 
56 
(56) 
cx = lim. ( > (c21)m) 
m1 


FINAL REMARKS 


By large x-values we should expect a zone of quickly damped disturbances in the 
neighborough of the origin & = 1 = 0. The Nx, Ny forces on the inner boundaries 
outside this zone might go nearly undisturbed to the outer boundaries with corres- 
pondingly bad convergence of the iteration procedure. It should be mentioned here 
that the nearly undamped solution may be somewhat incorrect, as the Donell solution 
for the boundary disturbances from the curved edge in cylindrical shells. ‘This is shown 
by Holand [5]. Hence, the error in the solution by repeated reflection of bondary 
disturbances may be considerable. However, the Nx and Ny forces on the outer bound- 
aries may be changed into shearing forces and moments by proper addition of solutions 
as the following 


U = K [ig — (57) 


where K is a real constant. From this equation and eq (11) we obtain with » = 0 
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Mx = 4K Dé# (58) 
M, =0 

= 16K Dg 

Qy =0 


The Nx, Ny forces which have to be eliminated with solutions of the same kind as 
eq (57), originate from the load of the inner stiffening beams and may for small beams 
be negligible. However, for a shallow shell loaded with a large snow load we need for 
longer spans a large edge beam to transmit the compressive forces. ‘The cross section 
can par example be made as shown in fig. 4. If the stiffness of the edge beam in the 
transverse direction is large compared to the stiffness of the shell and the width is 
small compared to the length, we can still use the boundary conditions eq (38). 

Considering now the shells with a x-value less than 2. For a shell with a = b and 
vy = 0 we obtain from eq (19) 


< 4,62 (59) 
Fig. 4. 


which corresponds to very shallow shells with comparatively short spans. However, 
if hollow sections (fig. 5) are used, shells with x < 2 covering great areas may be con- 
structed. One may realize this by rewriting eq (19) as 


yi— (60) 


where r is the radius of gyration of the cross section. 
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